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Abstract 



We report calculations of energy levels, radiative rates and electron impact excitation cross sections and 
rates for transitions in He-like Kr XXXV. The GRASP (general-purpose relativistic atomic structure package) 
is adopted for calculating energy levels and radiative rates. For determining the collision strengths and 
subsequently the excitation rates, the Dirac Atomic R-matrix Code (darc) is used. Oscillator strengths, 
radiative rates and line strengths are reported for all El, E2, Ml, and M2 transitions among the lowest 
49 levels. Additionally, theoretical lifetimes are listed for all 49 levels. Collision strengths are averaged 
over a Maxwellian velocity distribution and the effective collision strengths obtained listed over a wide 
temperature range up to 10 8 ' 1 K. Comparisons are made with similar data obtained with the Flexible 
Atomic Code (fac) to assess the accuracy of the results and to highlight the importance of resonances, 
included in calculations with DARC, in the determination of effective collision strengths. Differences between 
the collision strengths from DARC and FAC, particularly for forbidden transitions, are also discussed. Finally, 
discrepancies between the present results of effective collision strengths from the DARC code and earlier 
semi-relativistic .R-matrix data are noted over a wide range of electron temperatures for many transitions 
of Kr XXXV. 
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1 Introduction 



Krypton is used in a variety of experiments in tokamak fusion plasmas and has diagnostic applications. High 
temperatures in fusion plasmas, in excess of 10 7 K, give rise to many ionisation stages, including He-like Kr 
XXXV. Since krypton is increasingly being injected in fusion machines for the diagnostics of plasmas p], 
its study has become more important, particularly with the upcoming ITER project. However, to analyse 
observations, atomic data are required for a variety of parameters, such as energy levels, radiative rates (A- 
values), and excitation rates or equivalently the effective collision strengths (T), which are obtained from the 
electron impact collision strengths (f2). With this in view we have already reported atomic data for B-like and 
F-like Kr ions [2], [3], [4], and here we provide similar data for He- like Kr XXXV. 

Emission lines of He-like ions have been widely observed from a variety of laboratory plasmas. For example, 
spectra of Kr XXXV have been recorded in a high temperature Z-pinch discharge by Golt's et al [5]. Of 
particular interest are the resonance (w. Is 2 1 So - ls2p 1 P°), intercombination (x and y: Is 2 1 Sq - ls2p 
3 P2 ! i), and forbidden (z: Is 2 - ls2s 3 Si) lines of Kr XXXV [6], [7]. Considering the importance of Kr 
XXXV, a few calculations for its atomic parameters have already been performed using a variety of methods. 
Saloman [5] has compiled data from multiple calculations and has critically evaluated energy levels for many Kr 
ions, including Kr XXXV, and has posted data on the NIST (National Institute of Standards and Technology) 
website [http : //www . nistTg ov/pml/data/as d . cf mj Similarly, wavelengths are available for some transitions 
on the NIST website, but their quoted accuracy is not high. Recently, Zhang et al [9] have reported wavelengths 
and A- values, but only for the Is 2 1 So - 1 ' 3 P° transitions of Kr XXXV. Similarly, the collisional atomic data 
for Kr XXXV are also very limited. For example, Pindzola and Carter [10] have calculated collision strengths 
using the relativistic distorted-wave (DW) method, but only for transitions within the levels of the n=2 
configuration, and within a limited energy range of ~ 13-15 keV. Sampson et al |llj have reported collision 
strengths, but only for transitions from the lowest three levels to higher excited levels. Their calculations are 
based on the Coulomb-Born-exchange method and do not include the contribution of resonances, which can be 
very important as demonstrated in our earlier papers on other Kr ions [3], [I] as well as other He-like ions - see, 
for example, Aggarwal and Keenan [12] and references therein. Finally, Griffin and Ballance [13] have recently 
performed fully relativistic i?-matrix calculations for all transitions among the lowest 49 fine-structure levels 
of Kr XXXV, which belong to the Is 2 , ls2£, ls3£, ls4J, and ls5£ configurations, but have reported results of 
effective collision strengths only for transitions from the ground level up to ls4f F3, i.e. level 31. Therefore, 
in this paper we report a complete set of results (namely energy levels, radiative rates, lifetimes, and effective 
collision strengths) for all 1176 transitions among the lowest 49 fine-structure levels of Kr XXXV. Furthermore, 
we also provide the A- values for four types of transitions, namely electric dipole (El), electric quadrupole 
(E2), magnetic dipole (Ml), and magnetic quadrupole (M2), as these are required in a complete plasma model. 

For our calculations we employ the fully relativistic GRASP (general-purpose relativistic atomic structure 
package) code for the determination of wavefunctions, originally developed by Grant et al [14] and subse- 
quently revised by several workers, under the names GRASPl [15], GRASP2 [16], and GRASP2K [17] , How- 
ever, the version adopted here is GRASpO, which is based on [14] and is revised by Dr P H Norrington. 
This version contains most of the modifications undertaken in the other revised codes and is available on 
the website |http : //web . am . qub . ac . uk/DARCT] GRASP is a fully relativistic code, and is based on the jj 
coupling scheme. Further relativistic corrections arising from the Breit interaction and QED effects (vac- 
uum polarization and Lamb shift) have also been included. Additionally, we have used the option of ex- 
tended average level (EAL), in which a weighted (proportional to 2j+l) trace of the Hamiltonian matrix 
is minimized. This produces a compromise set of orbitals describing closely lying states with moderate 
accuracy. For our calculations of 51, we have adopted the Dirac Atomic R-matrix Code (darc) of P H 
Norrington and I P Grant (http://web.am.qub.ac.uk/DARC/). Finally, for comparison purposes we have 
performed parallel calculations with the Flexible Atomic Code (fac) of Gu [18], available from the website 
http : //sprg . ssl . berkeley . edu/~mf gu/f ac/. This is also a fully relativistic code which provides a variety of 
atomic parameters, and (generally) yields results for energy levels and radiative rates comparable to GRASP - 
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see, for example, [_], (T2], []_5]. However, differences in collision strengths and subsequently in effective collision 
strengths with those obtained from DARC can be large, particularly for forbidden transitions, as demonstrated 
in some of our earlier papers [12], and also discussed below in sections 5 and 6. Hence results from FAC will be 
helpful in assessing the accuracy of our energy levels and radiative rates, and in estimating the contribution 
of resonances to the determination of effective collision strengths, included in calculations from DARC but 
neglected in FAC. 

2 Energy levels 

The Is 2 , ls2£, ls3£, lsA£, and ls5£ configurations of Kr XXXV give rise to the lowest 49 levels listed in Table 
1, where we compare our level energies from GRASP (obtained without and with the inclusion of Breit and QED 
effects) with the critically compiled data by Saloman [8]. Also included in the table are our results obtained 
with the FAC code (FAC1), including the same CI (configuration interaction) as in GRASP. Our level energies 
obtained without the Breit and QED effects (GRASP1) are consistently higher than the NIST values by ~1.7 
Ryd. The inclusion of Breit and QED effects (GRASP2) affects (lowers) the energies by a maximum of ~2.1 
Ryd, and thus clearly demonstrates the importance of higher relativistic effects for a heavy ion, such as Kr 
XXXV. In addition, the orderings have slightly altered in a few instances, see for example the 4/5, 13/14 and 
23/24 levels. However, the energy differences for these swapped levels are very small. In general, the orderings 
from our GRASP calculations are nearly the same as those of NIST, but energies obtained with the inclusion 
of the Breit and QED effects (GRASP2) are consistently lower than the NIST values by ^0.3 Ryd, similar to 
the effect observed for other He-like ions [12] . Our FAC1 level energies are consistently higher by ^0.16 Ryd 
than the GRASP results, and hence are comparatively in better agreement with the NIST listings. Differences 
in the GRASP and FAC energies arise mostly by the ways calculations of central potential for radial orbitals 
and recoupling schemes of angular parts have been performed. Nevertheless, the level orderings from FACI 
are also in general agreement with our calculations from GRASP, but differ in some instances, particularly 
for the n = 5 levels. This is mainly because the degeneracy among the levels of the n = 5 configurations is 
very small. A further inclusion of the ls6£ configurations, labelled FAC2 calculations in Table 1, makes no 
appreciable difference either in the magnitude or orderings of the levels, mainly because the levels of the ls6£ 
configurations lie above the lowest 49 levels listed in Table I, and hence do not strongly interact with those. 

Finally, in Table 1 we also list the unpublished energies of Whitcford et al , which are obtained from 
the AutoStructure (as) code of Badnell [__], and are available at the APAP (Atomic Processes for Astrophysical 
Plasmas) website: http://amdpp.phys.strath.ac.uk/UK_APAP/. The AS energies are higher by up to ~2 
Ryd than the corresponding reference values from NIST, and our theoretical results from the FAC and GRASP 
codes. This is due to some of the relativistic effects being neglected in the AS calculations. More importantly, 
the level orderings are slightly different particularly for levels of the n = 5 configurations, and the AS energies 
for some of the levels are non-degenerate. However, the energy differences among the degenerate levels of a 
configuration are very small, as noted above. Since the NIST energies are not available for some of the levels, 
particularly of the n = 5 configurations, our energy levels either from the GRASP2 or FACI calculations 
should be adopted in modelling applications. 

3 Radiative rates 

The absorption oscillator strength (fij) and radiative rate Aji (in s _1 ) for a transition i — > j are related by 
the following expression: 

h = ^4^* = L49 x io- 16 a^M)^ (i) 

where m and e are the electron mass and charge, respectively, c is the velocity of light, Xji is the transition 
energy /wavelength in A, and uii and uij are the statistical weights of the lower (i) and upper (j) levels, 
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respectively. Similarly, the oscillator strength £y (dimcnsionless) and the line strength S (in atomic unit, 1 
a.u. = 6.460xl0 -36 cm 2 esu 2 ) are related by the standard equations given below. 

For the electric dipole (El) transitions 

2.02.1 x 10" ^ and flj = mil s ™ (2) 

for the magnetic dipole (Ml) transitions 

2.6974 x 10 13 _ M1 , , 4.044 x 10" 3 
—3 S and f i:j = — 

for the electric quadrupole (E2) transitions 

1.1199 xlQ 18 cE2 167.89 cE2 

A '* = T7T5 s and fa = -^T— s ' ( 4 ) 



A 3t = —-- S M1 and ^= , " S M \ (3) 
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and for the magnetic quadrupole (M2) transitions 

M.10X10» SM , ^ 2™_ S - (5) 

In Table 2 we present transition energies/wavelengths (A, in A), radiative rates (Aj i: in s _1 ), oscillator 
strengths (£y, dimensionless), and line strengths (S, in a.u.), in length form only, for all 336 electric dipole 
(El) transitions among the 49 levels of Kr XXXV considered here. The indices used to represent the lower and 
upper levels of a transition have already been defined in Table 1. Similarly, there are 391 electric quadrupole 
(E2), 316 magnetic dipole (Ml), and 410 magnetic quadrupole (M2) transitions among the 49 levels. However, 
for these transitions only the A-values are listed in Table 2, and the corresponding results for f- or S- values 
can be easily obtained using Eqs. (1-5). 

To assess the accuracy of our A-values, we have performed another calculation with the FAC code of Gu [TS] , 
For all strong transitions (f > 0.01), the A-values from GRASP and FAC agree to better than 10%. Furthermore, 
for a majority of the strong El transitions (f > 0.01) the length and velocity forms in our GRASP calculations 
agree to within 10%. However, the differences are larger for a few transitions, which are among the degenerate 
levels of a configuration, such as 4-5, 27-29 and 45-46, because their energy (AE) is very small and hence 
a slight variation in AE has a considerable effect on the A-values. For such transitions the two forms of the 
f- value differ by up to two orders of magnitude. Finally, as for the energy levels the effect of additional 
CI is negligible on the A- values, as results for all transitions agree within 10% with those obtained with 
the additional inclusion of the n = 6 configurations. To conclude, we may state that for almost all strong 
El transitions, our radiative rates are accurate to better than 10%. However, for the weaker transitions the 
accuracy is comparatively poorer. 



4 Lifetimes 

The lifetime r for a level j is defined as follows: 
1 



(0) 



Since this is a measurable parameter, it provides a check on the accuracy of the calculations. Therefore, 
in Table 1 we have also listed our calculated lifetimes, which include the contributions from four types of 
transitions, i.e. El, E2, Ml, and M2. To our knowledge, no calculations or measurements are available for 
lifetimes in Kr XXXV with which to compare. However, we hope the present results will be useful for future 
comparisons and may encourage experimentalists to measure lifetimes, particularly for the level ls2s which 
has comparatively a larger value. 
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5 Collision strengths 



Collision strengths (ft) are related to the more commonly known parameter collision cross section (cr^, 7rao 2 ) 
by the following relationship: 

n ij {E)=k}u) i a ij {E) (7) 

where kf is the incident energy of the electron and oji is the statistical weight of the initial state. Results for 
collisional data are preferred in the form of ft because it is a symmetric and dimensionless quantity. 

For the computation of collision strengths ft, we have employed the Dirac atomic R-matrix code (darc), 
which includes the relativistic effects in a systematic way, in both the target description and the scattering 
model. It is based on the jj coupling scheme, and uses the Dirac-Coulomb Hamiltonian in the i?-matrix 
approach. The i?-matrix radius adopted for Kr XXXV is 2.08 au, and 60 continuum orbitals have been 
included for each channel angular momentum in the expansion of the wavefunction, allowing us to compute 
ft up to an energy of 4000 Ryd, sufficient to determine values of effective collision strengths T (see section 6) 
up to T e = 10 81 K. The maximum number of channels for a partial wave is 217, and the corresponding size 
of the Hamiltonian matrix is 13 076. To obtain convergence of ft for all transitions and at all energies, we 
have included all partial waves with angular momentum J < 40.5, although a larger number would have been 
preferable for the convergence of some allowed transitions, especially at higher energies. However, to account 
for higher neglected partial waves, we have included a top-up, based on the Coulomb-Bcthe approximation 
[2"2] for allowed transitions and geometric series for others. 

For illustration, in Figs. 1-3 we show the variation of ft with angular momentum J for three transitions of 
Kr XXXV, namely 2-5 (ls2s 3 Si - ls2p 3 PJ), 2-11 (ls2s 3 Si - ls3p 3 P° 1 ), and 9-12 (ls3p 3 Pg - ls3p 3 P^), and 
at five energies of 1300, 1800, 2300, 2800 and 3300 Ryd. Values of ft have nearly converged for all resonance 
transitions (including the allowed ones), and most of the allowed transitions among the higher excited levels, 
as shown in Fig. 2 for the 2-11 transition. It is also clear from Fig. 2 that the need to include a larger range 
of partial waves increases with increasing energy. However, values of ft have not converged for those allowed 
transitions whose AE is very small (mainly within the same n complex), as shown for the 2-5 transition in 
Fig. 1. Similarly, values of ft have (almost) converged for all forbidden transitions, including those whose AE 
is very small, such as the 9-12 transition shown in Fig. 3. Therefore, mainly for the allowed transitions within 
the same n complex, our wide range of partial waves is not sufficient for the convergence of ft, for which a 
top-up has been included as mentioned above, and has been found to be appreciable. 

In Table 3 we list our values of ft for resonance transitions of Kr XXXV at energies above thresholds. The 
indices used to represent the levels of a transition have already been defined in Table 1. Unfortunately, no 
similar data are available for comparison purposes as already stated in section 1. Therefore, in order to make 
an accuracy assessment for ft, we have performed another calculation using the FAC code of Gu [15] . This 
code is also fully relativistic, and is based on the well-known and widely-used distorted-wave (DW) method. 
Furthermore, the same CI is included in FAC as in the calculations from DARC. Therefore, also included in 
Table 3 for comparison purposes are the ft values from FAC at a single excited energy Ej , which corresponds to 
an incident energy of ~ 2800 Ryd. For a majority of transitions the two sets of ft generally agree well (within 
~ 20%). However, the differences are larger for a few (particularly weaker) transitions. For example, for 64% 
of the Kr XXXV transitions, the values of ft agree within 20% at an energy of ~ 2800 Ryd, and discrepancies 
for others are mostly within a factor of two, although for some transitions (such as: 19-49, 33-38/41/43/47/49 
and 34-36/46/48), the differences are up to an order of magnitude. However, most of these transitions are 
weak (ft < 10~ 6 ) and forbidden, i.e. the values of ft have fully converged at all energies within our adopted 
range of partial waves in the calculations from the DARC code. For such weak transitions, values of ft from 
the FAC code are not assessed to be accurate. Additionally, for a few transitions values of ft from FAC are 
anomalous, as also noted for other He- like ions and demonstrated in Fig. 6 of Aggarwal and Keenan [T2], [23] . 
The sudden anomalous behaviour in values of ft from FAC is also responsible for the differences noted above 
for many of the transitions. Such anomalies for some transitions (both allowed and forbidden) from the FAC 
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calculations arise primarily because of the interpolation and extrapolation techniques employed in the code, 
which is designed to generate a large amount of atomic data in a comparatively very short period of time, 
and without too large loss of accuracy. Similarly, some differences in fi are expected because the DW method 
generally overestimates the results due to the exclusion of channel coupling. 

As a further comparison between the DARC and FAC values of fi, in Fig. 4 we show the variation of with 
energy for three allowed transitions among the excited levels of Kr XXXV, namely 2-6 (ls2s 3 Si - ls2p 3 P2), 5- 
14 (ls2p 3 P° - ls3d 3 D 2 ), and 11-24 (ls3p 3 P° - ls4d 3 D 2 ). Also included in this hgure are the corresponding 
results obtained with the FAC code. For many transitions there are no discrepancies between the f- values 
obtained with the two different codes (GRASP and fac), and therefore the values of ft also agree to better than 
20%. However, the values of ft obtained with FAC differ from our calculations with DARC, particularly towards 
the lower end of the energy range. Similar comparisons between the two calculations with DARC and FAC are 
shown in Fig. 5 for three forbidden transitions of Kr XXXV, namely 2-8 (ls2s 3 Si - ls3s 3 Si), 2-16 (ls2s 
3 Si - ls3d 3 D3), and 5-11 (ls2p 3 P° - ls3p 3 P°). As in the case of allowed transitions, for these (and many 
other) forbidden transitions the agreement between the two calculations improves considerably with increasing 
energy but the differences are significant towards the lower end of the energy range. These anomalies are 
due to the interpolation and extrapolation techniques employed in the FAC code, as stated above. Therefore, 
on the basis of these and other comparisons discussed above, collision strengths from the FAC code are not 
assessed to be very accurate, over the entire energy range, for a majority of transitions of Kr XXXV. However, 
we do not see any apparent deficiency in our calculations for ft, and estimate our results to be accurate to 
better than 20% for a majority of the (strong) transitions. 



6 Effective collision strengths 

Excitation rates, in addition to energy levels and radiative rates, are required for plasma modelling, and are 
determined from the collision strengths (ft). Since the threshold energy region is dominated by numerous 
closed-channel (Feshbach) resonances, values of ft need to be calculated in a fine energy mesh in order to 
accurately account for their contribution. Furthermore, in a plasma electrons have a wide distribution of 
velocities, and therefore values of ft are generally averaged over a Maxwellian distribution as follows: 

/>oo 

T(T e )= / fl(E)exp(-E j /kT e )d(E j /kT e ), (8) 
Jo 

where k is Boltzmann constant, T e is electron temperature in K, and E,,- is the electron energy with respect to 
the final (excited) state. Once the value of T is known the corresponding results for the excitation q(i,j) and 
de-excitation q(j,i) rates can be easily obtained from the following equations: 

q(i,j) = — Z Texp(-^/fcT e ) cmV 1 (9) 

UJiT e ' 



and 

, , 8.63 xlO" 6 ^ , , 

ffCj.O = ^T72— T cms > ( 10 ) 



0JjT e 

where Wi and tjj are the statistical weights of the initial (i) and final (J) states, respectively, and E,-j is the 
transition energy. The contribution of resonances may enhance the values of T over those of the background 
collision strengths (CIb), especially for the forbidden transitions, by up to an order of magnitude (or even more) 
depending on the transition and/or the temperature. Similarly, values of ft need to be calculated over a wide 
energy range (above thresholds) in order to obtain convergence of the integral in Eq. (8), as demonstrated in 
Fig. 7 of Aggarwal and Keenan [24] . 

To properly delineate resonances, we have performed our calculations of ft at over ~ 135 000 energies in 
the thresholds region. Close to thresholds (~0.1 Ryd above a threshold) the energy mesh is 0.001 Ryd, and 
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away from thresholds is 0.002 Ryd. Thus care has been taken to include as many resonances as possible, and 
with as fine a resolution as is computationally feasible. The density and importance of resonances can be 
appreciated from Figs. 7-10, where we plot 17 as a function of energy in the thresholds region for the four 
most important transitions of Kr XXXV, namely 1-2 (z: Is 2 1 S - ls2s 3 Si), 1-5 (y: Is 2 % - ls2p 3 P?), 1-6 
(x: Is 2 1 S Q - ls2p 3 P£), and 1-7 (w: Is 2 % - ls2p l PJ). For some transitions, such as 1-2, 1-5 and 1-6, 
the resonances are dense, particularly at energies just above the thresholds. These near threshold resonances 
affect the values of T particularly towards the lower end of the temperature range. Similar dense resonances 
have been noted for transitions in other He- like ions [12] . |23j . 

Our calculated values of T are listed in Table 4 over a wide temperature range up to 10 81 K, suitable 
for applications in a variety of plasmas. As stated in section 1, there are only limited results available for 
comparison purposes. Zhang and Sampson [25] have reported values of T for transitions among the lowest 
seven levels of Kr XXXV. In their calculations, they have adopted the Coulomb-Born-exchange method and 
have also included the contribution of resonances in an approximate way - see |25j for details and discussion. 
Therefore, in Table 5 we compare our results of T with theirs at three common temperatures of 1.87xl0 7 , 
4.20 xlO 7 and 9.33 xlO 7 K. Generally, the agreement between the two sets of results is within about a factor 
of two, such as for 4-5 (ls2p 3 P^ - ls2s 1 So). This is an 'elastic' transition (AE < 0.1 Ryd), has a very small 
f- value (^10 -5 ), and therefore the results of f2 can be highly variable as also noted for other ions. 

Since the above comparisons are limited to only a few transitions of Kr XXXV, we have also calculated 
values of T from our non-resonant f2 data obtained with the FAC code. These calculations are particularly 
helpful in providing an estimate of the importance of resonances in the determination of excitation rates. In 
Table 4 we have included these results from FAC at the lowest and the highest calculated temperatures for 
Kr XXXV. At T e = 10 67 K, our resonances-resolved values of T are higher by over 20% for about 45% of 
the transitions of Kr XXXV. Generally, the differences for a majority of the transitions are within a factor 
of two, but are higher (up to an order of magnitude) for a few, such as 18-28/29/30/31, 20-29/31, 24-30 
and 25-26. Most of these transitions are forbidden, and the differences are partly due to the contribution of 
resonances included in our calculations with the darc code. In the case of the most important w, x, y, and z 
lines, resonances have enhanced the values of T by about a factor of three for the z (1-2: Is 2 1 So - ls2s 3 Si) 
transition. A similar comparison at the highest temperature of our calculations, i.e. 10 81 K, indicates that 
about 40% of the transitions of Kr XXXV show differences of over 20% between the DARC and FAC values for 
T. These differences are generally within a factor of two, but are higher (up to a factor of eight) for a few, 
such as: 37-45 (ls5d 3 Di - ls5f 1 F%), 37-47 (ls5d 3 Di - ls5g 3 G 4 ), 38-44 (ls5d 3 D 2 - ls5f 3 F|) and 38-48 
(ls5d 3 D2 - ls5g 3 G5). All these four (and many other) transitions are forbidden, and the T from FAC are 
higher, mainly because of the sudden anomalies in the values of f2, as discussed in section 5. 

Based on the comparisons made above with the limited published results for transitions in Kr XXXV and 
with our FAC calculations, we estimate that the accuracy of our values of T from the DARC code listed in Table 
4 is better than 20% for a majority of transitions. A similar conclusion was also made for other He-like ions 
|12j . However, note that we have not included the effect of radiation damping in our calculations, which may 
reduce the importance of resonances in the collision strengths, for some of the transitions, and subsequently in 
the determination of T. This effect is important for highly charged ions, because radiative decay rates are large 
and compete with autoionisation rates [39] . However, earlier R- matrix calculations performed by Delahaye et 
al [39 with the inclusion of radiation damping show that their contribution (i.e. reduction in values of T) is 
appreciable only at lower temperatures (below 10 6 K), which may not be important for the plasma modelling 
for a highly ionised specie, such as Kr XXXV. Furthermore, Whiteford et al [20^ have demonstrated that the 
resultant uncertainties in the parts of the effective emission coefficients driven by excitation from the ground 
level for the four important lines of Ar XVII and Fe XXV (i.e. w, x, y and z), at appropriate temperatures 
and densities, is under 10%. Therefore, it is reasonable to conclude that our results of T listed in Table 4 
will not be significantly affected by the inclusion of radiation damping, and can be confidently applied to the 
modelling of plasmas. This in fact has been confirmed by the recent calculations of Griffin and Ballance |13) . 
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who performed similar calculations with the DARC code, with and without the inclusion of radiation damping, 
and concluded that the difference between the damped and undamped T values averaged over six temperatures 
(1.25 xlO 6 to 3.12xl0 8 K) for 30 resonance transitions is only 4.4%, and when averaged over 9 temperatures 
in the range 3.0xl0 6 to l.OxlO 9 K for all 1176 transitions of Kr XXXV, it is only 1.3%. 

In Table 6 we compare our results of T with the corresponding radiation damped values of Griffin and 
Ballance j!3j for the lowest 30 resonance transitions, which are common between the two calculations. In 
general, there is a good agreement between the two independent but similar calculations, although differences 
are up to ~20% for some of the transitions, particularly towards the lower end of the temperature range, 
which is not very important from the applications point of view. For some transitions, such as 1-8/13/14/17, 
our results are higher, whereas for some others, such as 1-28/29/30/31, those of Griffin and Ballance are 
higher. Such small differences between any two calculations are not uncommon and mainly arise because of 
the corresponding differences in the ways the calculations have been performed. For example, Griffin and 
Ballance calculated ft values up to an energy of 2000 Ryd and then extrapolated the data to higher energies 
in order to calculate corresponding values of T, whereas our calculations for fl are up to E = 4000 Ryd. These 
differences particularly affect the values of T at higher temperatures, and can be better appreciated in Fig. 
12 of Aggarwal and Keenan P2j and [25] . Similarly, Griffin and Ballance divided their calculations in three 
parts for three different ranges of partial waves, i.e. J <11.5, 12.5 < J < 30.5, and J > 31.5, whereas we 
performed a single calculation for all partial waves. Their procedure sometimes leads to discontinuity at the 
joining points for some of the transitions, and particularly affect the forbidden transitions. Finally, their energy 
mesh was comparatively coarser, .i.e. 0.00673 Ryd in comparison to ours of better than 0.002 Ryd, and they 
had even a broader spacing of 0.648 Ryd for partial waves with J > 12.5, whereas we have a uniform energy 
mesh for all partial waves. These differences in the energy spacings particularly affect the determination of the 
contribution of resonances. Considering these differences in the calculation procedures, differences in T values 
noted in Table 6 are fully understandable. However, these comparisons are limited to only ^2.5% transitions 
and therefore we discuss below the comparisons among a wider range of transitions. 

As for the energy levels listed in Table 1, effective collision strengths for transitions in Kr XXXV are 
available on the APAP website: http : //amdpp . phys . strath . ac . uk/UK_APAP/. These results were undertaken 
alongside the earlier calculations of Whiteford et al [20] for Ar XVII and Fe XXV. For the calculations, they 
adopted the standard i?-matrix code of Berrington et al |40j . The electron exchange was included for partial 
waves up to J = 10.5 but was neglected for higher J values. Furthermore, the calculations were performed in 
the LS coupling scheme and the corresponding results for 51 and subsequently T for fine-structure transitions 
were obtained using their intermediate coupling frame transformation (icft) method. However, the data 
obtained by such procedures should qualitatively be comparable with our fully relativistic results from the 
DARC code, as has already been demonstrated in several papers - see for example, Liang and Badnell [41] . 
Furthermore, Whiteford et al included the effect of radiation damping, which can reduce the contribution of 
resonances in the determination of T values, for some transitions, such as ls2p 1 PJ - ls3s ^o, particularly 
towards the lower end of the temperature range - see their Fig. 4. However, at temperatures relevant to 
plasma modelling applications for Kr XXXV, the effect of radiation damping is negligible as discussed above. 

As stated above, the two independent sets of results from DARC and Whiteford et al [20\ should be broadly 
comparable, although relativistic effects are very important for a heavy ion considered in this paper, and 
noted earlier for energy levels in section 2. However, differences between the two datasets are considerable 
(over 20%) for many transitions, and throughout the temperature range of the calculations. For example, 
at the lowest common temperature of 2.45 xlO 5 K, the two sets of T differ by over 20% for ~40% of the 
transitions. For a majority of transitions, these differences are within a factor of three. For some transitions 
our results are higher, but for most others the reverse is true. Of particular note are those transitions for 
which the differences are over two orders of magnitude, and the T values of Whiteford et al are invariably 
higher. Examples of such transitions are: 13-14/16/17, 14-16/17 and 16-17. Most of these transitions belong 
to the degenerate levels of a state/configuration, and hence have very small transition energies, as seen in Table 
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1. Such transitions may be referred to as 'elastic' and are always problematic to calculate - see for example 
Hamada et al [32] for hydrogenic ions which show similar problems. To demonstrate the differences between 
the two calculations, in Fig. 10 we compare values of T for three transitions of Kr XXXV, namely 13-14 (ls3d 
3 Di - ls3d 3 D 2 ), 14-16 (ls3d 3 D 2 - ls3d 3 D 3 ), and 29-31 (ls4f 3 F§ - ls4f 1 F%). These differences in the T 
values are not due to any resonances, but arise from the limitation of the approach adopted by Whiteford et al 
|20j . as recently discussed and demonstrated by Bautista et al. [43]. Similar large differences are also observed 
with the calculations of Whiteford et al [H] for transitions of Li-like ions, as discussed and demonstrated by 
Aggarwal and Keenan [45], [46]. The problem in the i?-matrix code adopted by Whiteford et al [20], [44] has 
been identified and rectified by Liang and Badnell |41j . However, since the calculations of Whiteford et al [20] 
were performed more than a decade ago, limitations in their data for T remain, and a reexamination of their 
results is desirable. 

The differences between our results of T from DARC and those of Whiteford et al [20] are not confined 
to lower temperatures, but are also found at higher temperatures. At the highest common temperature of 
1.22xl0 8 K, the two sets of T differ by over 20% for -15% of the transitions of Kr XXXV. Hence, there 
is comparatively a better convergence of the results at higher temperatures. However, for almost all such 
transitions, the T values of Whiteford et al are higher. The differences for most transitions are within a factor 
of two, but for some are up to an order of magnitude, such as: 14-16/17, 19-43/45, 29-31 and 33-43/45. 
All these transitions are forbidden and Q for these have converged within our adopted partial waves range, as 
discussed in section 5. To demonstrate the differences, in Fig. 11 we compare our results of T from DARC 
with those of Whiteford et al [20] for three transitions of Kr XXXV, namely 14-17 (ls3d 3 Di - ls3d ^2), 
46-48 (ls5g 3 G 3 - ls5g 3 G 5 ) and 47-49 (ls5g 3 G 4 - ls5g 1 G 4 ). Since the collision strengths of Whiteford et 
al are overestimated for many transitions, the corresponding results for T are affected throughout the entire 
temperature range. 

To conclude, we may state that the T results of Whiteford et al [20 should not be as accurate as those 
presented here for the reasons given above. Based on several comparisons as well as with our experience with 
other He- like ions [12] , [23] and [33] , we estimate that the accuracy of our T results from the DARC code listed 
in Table 4 to be better than 20% for a majority of transitions. 

7 Conclusions 

In this paper we have presented results for energy levels and radiative rates for four types of transitions (El, 
E2, Ml, and M2) among the lowest 49 levels of Kr XXXV belonging to the n < 5 configurations. Additionally, 
lifetimes of all the calculated levels have been reported, although no measurements or other theoretical results 
are available for comparison. However, based on a variety of comparisons among various calculations with 
the GRASP and FAC codes, our results for radiative rates, oscillator strengths, line strengths, and lifetimes 
are judged to be accurate to better than 20% for a majority of the strong transitions (levels). Additionally, 
we have considered a large number of partial waves to achieve convergence of at all energies, included a 
wide energy range to accurately calculate values of T up to T e = 10 81 K, and resolved resonances in a fine 
energy mesh to account for their contributions. Hence we see no apparent deficiency in our reported results for 
collision strengths and effective collision strengths, and estimate their accuracy to be better than 20% for most 
transitions. However, the present data for effective collision strengths for transitions involving the levels of the 
n = 5 configurations may be improved somewhat by the inclusion of the levels of the n = 6 configurations. 
We believe the present set of complete results for radiative and excitation rates for all transitions of Kr XXXV 
will be highly useful for the modelling of a variety of plasmas. 
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Figure 1. Partial collision strengths for the 2 — 5 transition of Kr XXXV. 
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Figure 1: Partial collision strengths for the ls2s 3 Si - ls2p 3 P^ (2-5) transition of Kr XXXV, at five energies 
of: 1300 Ryd (circles), 1800 Ryd (triangles), 2300 Ryd (stars), 2800 Ryd (squares), and 3300 Ryd (diamonds). 
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Figure 2. Partial collision strengths for the 2—11 transition of Kr XXXV. 
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Figure 2: Partial collision strengths for the ls2s 3 Si - ls3p 3 P ! j > (2-11) transition of Kr XXXV, at five energies 
of: 1300 Ryd (circles), 1800 Ryd (triangles), 2300 Ryd (stars), 2800 Ryd (squares), and 3300 Ryd (diamonds). 
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Figure 3. Partial collision strengths for the 9-12 transition of Kr XXXV. 
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Figure 3: Partial collision strengths for the ls3p 3 Pq - ls3p (9-12) transition of Kr XXXV, at five energies 
of: 1300 Ryd (circles), 1800 Ryd (triangles), 2300 Ryd (stars), 2800 Ryd (squares), and 3300 Ryd (diamonds). 
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Figure 4. Comparison of collision strengths for the 2 — 6, 5—14, and 11—24 allowed transitions of Kr XXXV. 
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Figure 4: Comparison of collision strengths from our calculations from darc (continuous curves) and FAC 
(broken curves) for the 2-6 (circles: ls2s 3 Si - ls2p 3 P§), 5-14 (triangles : ls2p 3 P? - ls3d 3 D 2 ), and 11-24 
(stars: ls3p 3 P? - ls4d 3 D 2 ) allowed transitions of Kr XXXV. 
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Figure 5. Comparison of collision strengths for the 2—8, 2—16, and 5-11 forbidden transitions of Kr 
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Figure 5: Comparison of collision strengths from our calculations from darc (continuous curves) and FAC 
(broken curves) for the 2-8 (circles: ls2s 3 Si - ls3s 3 Si), 2-16 (triangles: ls2s 3 Si - ls3d 3 D 3 ), and 5-11 
(stars: ls2p 3 P? - ls3p 3 P1) forbidden transitions of Kr XXXV. 
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Fig. 6. Collision strengths for the 1 -2 transition of Kr XXXV 
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Figure 6: Collision strengths for the Is 2 % - ls2s 3 Si (1-2) transition of Kr XXXV. 
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Fig. 7. Collision strengths for the 1 -5 transition of Kr XXXV 

0.020 1 — n im r 



0.016 - 



0.012 



in 



c 
o 



•I 0.008 
o 
o 



0.004 - 



0.000 



950 



1000 



I Z 



1050 



1 100 
E (Ryd) 



T 



1 150 



T 



1200 



1250 



Figure 7: Collision strengths for the Is 2 % - ls2p 3 P? (1-5) transition of Kr XXXV. 



19 



Fig. 8. Collision strengths for the 1 -6 transition of Kr XXXV 
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Figure 8: Collision strengths for the Is 2 % - ls2p 3 P% (1-6) transition of Kr XXXV. 
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Figure 9: Collision strengths for the Is 2 % - ls2p ^ (1-7) transition of Kr XXXV. 
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Figure 10. Comparison of effective collision strengths for the 13—14, 14—16 and 29 — 31 transitions of Kr XXXV. 
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Figure 10: Comparison of effective collision strengths for the 13-14 (circles: ls3d 3 Di - ls3d 3 D2), 14-16 
(triangles: ls3d 3 D 2 - ls3d 3 D 3 ), and 29-31 (stars: ls4f 3 F° - ls4f 1 F^) transitions of Kr XXXV. Continuous 
and dotted curves are from the present DARC and earlier R- matrix codes [20], respectively. 
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Figure 11. Comparison of effective collision strengths for the 14—17, 46 — 48 and 47 — 49 transitions of Kr XXXV. 
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Figure 11: Comparison of effective collision strengths for the 14-17 (circles: ls3d 3 D2 - ls3d 1 T>2), 46-48 (ls5g 
3 G 3 - ls5g 3 G 5 ), and 47-49 (stars: ls5g 3 G 4 - ls5g 1 G 4 ) transitions of Kr XXXV. Continuous and dotted 
curves are from the present DARC and earlier R- matrix codes [20], respectively. 
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Table 1: Energy levels (in Ryd) of Kr XXXV and their lifetimes (r, s). a±b = axlO ±fc . 
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Tabic 3: Collision strengths for transitions in Kr XXXV. (a±b = axlO ±b ). 
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Table 5: Comparison of T values for transitions of Kr XXXV. (a±b = axlO ±h ). 
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DARC: Present calculations from the DARC code 
ZS: Calculations of Zhang and Sampson [25] 
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Table 6: Comparison of T values for resonance transitions of Kr XXXV. (a±b = axlO ±b ). 
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2.20- 


-6 


6.055- 


-7 


6.96 


-7 


3.707 


-7 


4.18- 


-7 


1 
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3.531 


-6 


4.32 


-6 


2.264- 


-6 


2.74- 


-6 


7.920 


-7 


9.04 


-7 


5.378- 


-7 


6.02 


-7 


1 


30 


3.743 


-6 


4.52 


-6 


2.454- 


-6 


2.93 


-6 


8.693- 


-7 


9.83 


-7 


5.422- 


-7 


6.02 


-7 


1 


31 


3.423 


-6 


4.35 


-6 


2.168 


-6 


2.73 


-6 


7.634 


-7 


8.94 


-7 


5.388- 


-7 


6.10- 


-7 



AK: Present calculations from the DARC code without radiation damping 

GB: Calculations of Griffin and Ballance [13] from the DARC code with radiation damping 
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